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' The specific heat coefficient "f s (H) and the electronic thermal conductivity n ea {H) are calculated 

for Abrikosov's vortex lattice by taking into account the effects of supercurrent flow and Andreev 
scattering. First we solve the gap equation for the entire range of magnetic fields. We take into 
account vertex corrections due to impurity scattering calculated in the Born approximation. The 
function j a (H)/y n increases from zero and becomes approximately linear above H/H C 2 ~ O.f . The 
dependence on impurity scattering is substantially reduced by the vertex corrections. The upward 
curvature of K es (H)/ K en , which is caused by decreasing Andreev scattering for increasing field, is 
reduced for increasing impurity scattering. We also calculate the temperature dependence of the 
scattering rates l/r ps (H) of a phonon and l/T ea (H) of a quasiparticle due to quasiparticle and 
phonon scattering, respectively. At low temperatures the ratio T pn /T ps (H) increases rapidly to one 
as H tends to H C 2 which yields a rapid drop in the phononic thermal conductivity K p h- Our results 
^ \ are in qualitative agreement with the experiments on the two-gap superconductor MgB2 . 
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£ ■ I. INTRODUCTION 

o 

Although superconductivity in MgB2 1 has been established as conventional phonon-mediated s-wave superconduc- 
tivity, it is often difficult to analyze the results of measurements of various physical quantities because there exist two 
gaps of different magnitude associated with different bands. Their ratio is estimated as r — Af /Aq - 0.3-0.4 where 
the larger gap Aq is associated with the two-dimensional cr-bands and the smaller gap Aq with the three-dimensional 
7r-bands. Evidence for two gaps is provided by the rapid rise of the specific heat coefficient J S {H) at very low fields. 2 ' 3 
' Recent measurements of the in-plane thermal conductivity n s (H) for fields both parallel and perpendicular to the 
c-axis provide more evidence for the existence of two gaps because they show a very unusual field dependence at low 
t-H ' temperatures. For increasing field, k s (H) drops rapidly and reaches a minimum at a relatively low field, then exhibits 
an S-shape behavior where a steep rise in the low-field region is followed by a downward curvature up to about H c2 /2 
and finally an upward curvature near H C 2 for H || c. The initial drop is attributed to the decrease of the phononic 
contribution K p h(H) which is caused by the rapid enhancement of scattering of phonons by the quasiparticles in the 
vortex cores associated with the 7r-band. The steep increase in the low field region is attributed to the rapid rise of 
the electronic contribution n e (H) which arises from the the rapid release of mobile quasiparticles in the vortex lattice 
. asociated with the 7r-band. The final upward curvature of k(H) near H c2 for H || c is presumed to be due to the 
c~j ' mobile carriers of the vortex lattice associated with the cr-band. 

The field dependence of "f s {H) and K es {H) can be explained qualitatively in terms of a two-band model with 
O ■ two energy gaps of different magnitude in the two bands where the band with strong pairing (the L, or cr-band) is 
responsible for the superconductivity and the superconductivity in the second band (the S, or 7r-band) is induced by 
Cooper pair tunneling. 5 The vortices in the S-band have large vortex cores of radius £q and already start to overlap 
in weak fields. This leads to a strong suppression of the small gap at a "virtual" upper critical field H^ 2 ~ 3 — 4 kOe 

' (c) 

which is much smaller than H^ 2 ~ 35 kOe. The large radius of the vortex core and the field-induced suppression of 
the smaller gap is consistent with recent scanning tunneling spectroscopy (STS) measurements. 6 

Theories of the electronic thermal conductivity n e 7 and the phononic thermal conductivity n p h S in the mixed 
state have been developed on the basis of the Brandt-Pesch- Tewordt (BPT) 9-11 and the equivalent Pesch (P) 12 
approximation schemes. Linear response equations for the vortex state 13 yield another expression for n e which 
has been applied to d-wave pairing superconductivity in the cuprates. 14 In the BPT-approximation the normal and 
anomalous Green's functions G and F are derived from Gorkov's integral equations with kernel A(ri)C?^_ w (ri — 

r2)A*(r2) exp(— i 2eA ■ ds) where A(r) is Abrikosov's vortex lattice order parameter and G°_ u is the quasihole 
propagator. Thus these Green's functions take into account the effects of supercurrent flow and Andreev scattering 
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in the vortex cores as well as quasiparticle transfer between the vortices. The expression for K e 7 , which has been 
derived from the Kubo formula in analogy to the procedure in Ref. 15, consists essentially of the w-integral with 
the factor (o;/T) 2 sech 2 (o;/2T) determining the temperature dependence, the "destructive" coherence factor arising 
from the combination (GG^ — FF^ ) , and the total scattering rate in the denominator consisting of the sum of the 
quasiparticle scattering rates 7 and 7,4(0;, i?) due to impurity and Andreev scattering. In the present paper we shall 
also calculate the scattering rate l/r es due to thermal scattering in the vortex state, but this becomes comparable to 
7 only at higher temperatures. The expression for K p h contains, in addition to the sum of scattering rates of phonons 
by different defects, the scattering rate l/r ps for phonons scattered by quasiparticles in the vortex state. The latter 
expression was obtained by replacing in the BCS-expression 16 the density of states (DOS) and "coherence" functions 
by the expressions resulting from the P-approximation for the vortex state. 8 

The results for -f s (H)/j n , n es {H)/n en , and r p „/r ps , which were obtained in Ref. 8 for a single cylindrical or 
spherical band and an isotropic s-wave pairing gap, turned out to depend sensitively on the reduced impurity scattering 
rate 5 = T/Ao = (it/2)£ > o/£ where £0 is the zero-field coherence length and £ is the mean free path. If one takes for 
the 7r-band in MgB2 £„■ ~ 500A from the STS measurements 6 and a mean free path £ n ~ 500 — 800A , one arrives at a 
reduced sacttering rate 8 V ~ 1 which is much larger than the largest value of 5 = 0.5 taken in Ref. 8. For such a large 
impurity scattering rate it is necessary to renormalize the impurity scattering self energy and to take into account the 
corresponding vertex corrections. The vertex corrections yield a renormalization of the gap function 10,12 which, near 
H C 2 , tends to the renormalization function which has been derived in the calculation of H C 2 . 17 It turns out that this 
gap renormalization substantially reduces the dependence of the functions j s (H)/'y n , K es(H) / K en , and T pn /T ps (H) 
on impurity scattering. 

Recently it has been shown that numerical solutions of the Eilenberger equations yield a spatial average of the 
density of states (DOS) which is approximated very well by the DOS obtained with the BPT- and P-approximations 
over the whole field range between H C 2 and H c \ . ls This result motivates us to calculate the field dependence of 
the central parameter of the theory, A = AA/u, from the gap equation in the P-approximation 18 by including the 
impurity scattering vertex corrections in analogy to those in the gap equation of the BPT- approximation. 10 Here A 2 
is the spatial average of |A(r)| 2 where A(r) is Abrikosov's vortex lattice order parameter, A = (2eH)~ 1 / 2 = (7r/2) 1 / 2 a 
where a is the lattice constant of the vortex lattice, and v is the Fermi velocity perpendicular to the uniform field 
H. It turns out that one can fit the exact results of these calculations very well with the relation between A 2 and 
H/H C 2 given by the Ginzburg-Landau relationship if one uses an appropriate value of the Abrikosov parameter 
(3 A =< |A| 4 > /(< |A| 2 >) 2 . 

We have also investigated the temperature dependence of the phonon relaxation time ratio T pn / r ps (H) and find that 
this ratio increases rapidly to one with increasing H in proportion to [y s (H)/^ n ] 2 for temperatures T < (1/2) A(T, H) . 
This yields a rapid decrease of n p h for small fields and relatively low temperatures if one takes for A the small gap 
A S (T,H) of the 7r-band. In this way one can explain the observed drop of k for small fields and temperatures near 
T c /6 which is attributed to the drop of K p h - 4 

Since the interband impurity scattering between the a- and 7r-band is small, 20 one can add, in a crude approximation, 
our results obtained for a single isotropic gap in a cylindrical or spherical Fermi surface (FS) band by weighting them 
with the corresponding density of states. If one takes into account that the "virtual" upper critical field for the three- 
dimensional 7r-band is much smaller than the upper critical field of the two-dimensional cr-band, 3 one can explain 
qualitatively the measured field dependence of the specific heat and thermal conductivity in MgB 2 . 



II. THEORY OF THE SPECIFIC HEAT AND ELECTRONIC AND PHONONIC THERMAL 
CONDUCTIVITY IN ABRIKOSOV'S VORTEX LATTICE 



The basis of our theory is the approximation of Brandt, Pesch, and Tewordt (BPT) 9-11 for the normal and anoma- 
lous Green's functions G and F in the mixed state. A simplified version of this theory has been derived by Pesch (P) 12 
from the quasiclassical Eilenberger equations. In this version the e-integrals of the spatial averages of the spectral 
functions, — ImG/n and — ImF/7r , denoted by A and B , take on the following form: 8 
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z = 2[QA + i(A/v)j]/sm9; 0=Z(p,H), 
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A = (2eH)~ 1/2 , A = AA/u , = w/A , (4) 
A 2 = (ff c2 - ff)/6/? A ff , A 2 = A 2 (T) [1 - (H/ffca)] , (5) 
w(z) = cxp(— z 2 )erfc(— iz) . (6) 

A r = A/ D , D(fl) = l- d9 2{A 1 /v)^w(z). (7) 
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Here A (T) is the s-wave gap in zero field, ff is the spatial average of the magnetic field, A 2 =< |A| 2 > , where A(r) 
is Abrikosov's vortex lattice order parameter, and Pa is the Abrikosov parameter. For the total impurity scattering 
rate 7 we employ the Born approximation, 7 = TA where T = l/2r„ and A is the average of A over the Fermi 
surface which is calculated self-consistently. For a spherical Fermi surface (FS), usin# is the component of the Fermi 
velocity v(p) perpendicular to the field H . For a cylindrical FS parallel to c and field H || c , one has to take 9 = n/2 
and v equal to the Fermi velocity in the ab-plane. For a general FS, vsinO is replaced by the component i>j_(p) 
perpendicular to the field. The gap renormalization function D in Eq.(7) was first calculated by U. Brandt 10 from 
the ladder summation of impurity interaction lines bridging the vertices at A and A* occuring in the self energy part 
G° AG°_ W A*G°_, in Gorkov's integral equation for G . In the P-approxiamtion scheme D takes the form given in Eq.(7) 
which, near H c2 , tends to the function introduced in Ref. 17. 

The field dependence of the DOS, A, and the "coherence factor" function, B, are governed by the quantity A = 
AA/v = {2n)~ 1 / 2 a/£, where a = {2/n) 1 / 2 A is the lattice constant of the Abrikosov vortex lattice and £ = v/Att is the 
effective coherence length in the plane perpendicular to H . For H — ► H C 2 and H —* one obtains the limiting values 
A — > and A — > 00 , respectively. The relationships between A and H/H C 2 and between A and H/H C 2 given in 
Eq.(5) have been derived from the GL-theory which is valid near H C 2 . To get the general relationship at lower fields 
we make use of the gap equation in the P-approximation. 18 At T = this equation can be rewritten in the following 
form: 

fOO 

dn [B(n) - = . (8) 

Here £?(Q) is the spectral function of F given in Eq. (2) and C(fi) is obtained by taking the limit of B for H — > H c2 
where A — > : 



C(fi) -Re -i^2(H/H c2 ) 1/2 (A r /sm8)w (H/H c2 ) 1/2 z . (9) 

L L J J A(H c2 ) 

The gap equation Eq. (8) in the P-approximation corresponds to the gap equation in the BPT-approximation 10 with 
the function B replaced by the spectral function of F derived for the theory of NMR. 11 Thus Eq. (8) is an implicit 
equation for H/H C 2 as a function of A for given impurity scattering rate 8 = T/Aq. This calculation is rather 
complicated because the impurity scattering self energy 7 = TA(fl) occuring in the expression for B(fl) in Eq. (2) 
has to be calculated self-consistently for each value of from the expression for A(Q) (sec Eq. (1)). For H || c and 
a cylindrical FS one has to set 9 = ir/2 in Eq.(8), and for a spherical FS Eq.(8) is averaged over the polar angle 9 . 
In Fig. I we show our results for H/H C 2 versus A obtained from Eq.(8) for 9 = n/2 and for 5 = 0.1, 0.5, and 1.0. 
These results can be fitted very well by the relationship H/H c2 = (1 + 6/3aA 2 ) -1 [see Eq.(5)] over the whole field 
range if one takes Pa — 1-1, 1.4, and 1.7 for 5 = 0.1, 0.5, and 1.0. Similar results are obtained for a spherical FS if 
one averages Eq.(8) over the polar angle 9 (see the dashed curves in Fig. 1 for S = 0.1, 0.5, and 0.8 with (3a — 2.3, 
2.8, and 3.2). 

In Fig. 2 we have plotted our numerical results for the normalized DOS at zero energy, A(Q, = 0) = -f s (H)/'j n , 
versus H/H C 2 where 7 is the specific heat coefficient. Here we have approximated A occuring in Eq. (1) by the 
analytical expression in Eq. (5) with (3a = 1-1, 1.4, and 1.7 for bare impurity scattering rates 8 — T / A = 0.1, 0.5, 
and 1.0 (solid curves, from top to bottom). Here we have employed the Born approximation for the total impurity 
scattering self energy, 7 = T A(0) , where A(0) is calculated self-consistently. This means that 8 is replaced by 8 A(0) . 
For 8 = 0.1 and 9 — n/2 the function A(0) is close to the clean limit result shown in Ref. 18 for the s-wave pairing 
state. The latter function has been shown to be very close to the numerical solution of the Eilcnbcrgcr equations over 



3 



the whole field range from H C 2 to H = corresponding to H c \ . One can see that all curves in Fig. 2 tend to zero 
in the limit H / H C 2 — > , which is due to the gap renormalization. For a spherical FS one has to average Eq.(l) over 
the polar angle 9 . For the plot of A versus H/H C 2 one needs Eq.(5) where Pa is calculated from the #-angle average 
of Eq.(8). The results are shown by the dashed curves in Fig. 2 for S = 0.1, 0.5, and 0.8 with the /^-values used in 
Fig. 1. It is interesting that the curve for S = 0.8 is almost linear over the whole field region. 

We turn now to the theory of the electronic thermal conductivity n e in the vortex state which has been developed 
in Ref. 7 by replacing the zero-field Green's functions G and F in the theory of Ref. 15 by the Green's functions of the 
BPT-approximation. This theory has been applied to superconducting states with nodes in the gap such as Sr 2 Ru04 , 
the cuprates, and UPt3 . To save space we omit here the expression for n e . The main features of this expression 
are the following. Most important is the term Imeo in the denominator of the w-integral with the well-known factor 
w 2 sech 2 (w/2T) . This term Ime corresponds to the scattering rate of quasiparticles due to impurity and Andreev 
scattering. The other term corresponds to the coherence factor of BCS theory. In the zero field limit both terms in 
the expression of Ref. 7 tend correctly to those occuring in the expression of Ref. 15. The physical meaning of Imeo 
is the following. The equation for the position, eo , of the pole of the Green's function G in the BPT-approximation 9 
yields Im eo =7 + 7.4 where 7 is the total impurity scattering rate and 7 a is the imaginary part of the quasiparticle 
self energy E w at eo • From the kernel of Gorkov's integral equation for G one sees that, in the spatial representation, 
£ w (ri,r 2 ) = — V(ri,r2)G°_ UJ (ri — r 2 ) , where V{ri,r 2 ) — A(ri)A*(r 2 ) exp(— 1 f^2eA ■ ds) and A(r) is Abrikosov's 
vortex lattice order parameter. From this expression one recognizes that ja = — Im£ w is the scattering rate for 
converting a quasiparticle into a quasihole by Andreev reflection at A*(r 2 ) and then back into a quasiparticle at 
A(n) . The Fourier transform of S w with respect to the difference coordinate (ri — r 2 ) is equal to 

(p) = - J d 3 p> V{ P ')G\ (p - p') , (10) 

with 

V(p') = 8n 2 A 2 A 2 5(p' z ) exp [-A 2 (p 12 + p' y 2 )} (z || H) . (11) 

For a sperical Fermi surface the result is 9 

= -i\A r A 2 (A/usin6»)ii;[(w + + e)A/v sm9] (12) 

where e is the normal state energy measured from the Fermi energy, 9 = Z(p, H) , and 7 is the total impurity scattering 
self energy. 

In the uo — > limit corresponding to T — > one obtains the following explicit expression for Imeo ;7 ' 10 
Ime = 7 + 1A ; 1A = iy/K)(3 

_ fr _ n 1/2 

13 = A 2 I [( 7 A/w) 2 + A 2 + 77 sin 2 9\ + (jA/v) \ ; (tt^ 1 < rj < 1/2) . (13) 

All the other quantities A and (3 occuring in the expression for K e (see Eq.(13) of Ref. 7)) are also renormalized by 
the replacement A — > A r (see Eq.(7) for fl = 0). Note that A <~ a/£ and TA/v ~ a/l . In Fig. 3 we have plotted our 
results for n es / n en versus H/ H C 2 in the limit T —t . The bare impurity scattering rates 5 = T/ Aq are 0.1, 0.5, and 
1.0 (from bottom to top). For a spherical FS one has to evaluate the average of n e over the polar angle 9 which yields 
nearly the same results as shown in Fig. 3 for the plots of K es /n en vs H/H C 2 if one uses Eq.(5) with the Pa calculated 
from the average of Eq.(8) (see the dashed curves in Fig. 3 for <5 = 0.1, 0.5, and 0.8 with the /3^-values given in Fig. 
1). All curves tend to zero in the limit H/H C 2 — » . One sees that, for low impurity scattering rates, n es / n en exhibits 
a strong upward curvature near iJ c2 which is caused by the rapid reduction of the Andreev scattering rate ja with 
decreasing A or increasing field as can be seen from Eq. (13). For increasing impurity scattering rate 7 this effect is 
diminished because ja becomes smaller in comparison to 7 as can be seen from the expression in Eq. (13). 

The main problem in analyzing the measured thermal conductivity k in the mixed state is to separate the electronic 
and phononic contributions n e and n p h . 4 The general expression for n p h (see Eq. (4) of Ref. 8) contains, in addition 
to other phonon relaxation times due to defect scattering, the relaxation time t p (uj) due to scattering of phonons of 
frequency uu by quasiparticles which changes drastically in the BCS state. The expression for l/r ps in the BCS state 16 
has been extended to the vortex lattice state by replacing the BCS DOS E/e , where e = (E 2 — A 2 ) 1 / 2 , and the 
function A/e occuring in the BCS coherence factor, with the spectral functions A and B of Eqs. (1) and (2). In the 
zero field limit the latter expressions tend correctly to the BCS functions with u> = E . The resulting ratio of phonon 
relaxation times in the normal and mixed states is given by (see Eq. (3) of Ref. 8): 
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T P n/T ps =fl(fi ) = [l-exp(-n /y)](2/O ) / dn/[(fi-fio/2)/tf]/[-(n + ft /2)/y] 
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Here, fl = wo/A is the phonon energy lu divided by A = A(T,H) , f is the Fermi function, and y — T/A(T,H) . 
In Ref. 8 we calculated g in the limit T — > . It was found that in the limit A — > oo , or H/ H C 2 — > , g tends 
correctly to the BCS step functon with a step at fl = 2 , or uj = 2A . This step function is washed out as A — ► , 
or H/H c2 — > 1 , and it tends to the constant, 1. Since it is important to know the temperature dependence of K p h , 
and thus of T pn /T ps — g , for analyzing the measured k as a function of T at fixed field, 4 we have now calculated 
the T-dependence of the expression in Eq.(14). For this purpose we make use of the fact that, according to the 
general expression for K p h (sec Eq.(4) of Ref. 8), the maximum of the integrand in that x-intcgral occurs at about 
x m = uJom/T ~ 4 . Therefore it is a good approximation to make the variable transformations Qq = x m y and Q = ty 
in Eq.(14). This yields t p „/t ps = g(x m y) as a function of the T-dependent quantity y for fixed A, or field H/H c2 . 
In Fig. 4 we have plotted g(x m y) versus y for constant x m — 4 and constant A = 0.2, 0.3, and 0.6, or H/H c2 = 0.78, 
0.62, and 0.29 (see Fig. 1) for 6 = 0.1 . In the limit H/H C 2 —* , one obtains a step function with a step at y = 1/2 , 
that is, f^o = Mj — 2 or uj = 2A , as it should be. One sees, that for intermediate values of the field, T pn /T ps = g 
considered as a function of T is smeared out and tends to one for high fields. In the limit y — > , or £Iq — > Eq.(14) 
yields g(Cl = 0) = L4(f2o = 0, A , 8)] 2 . If the scattering rate l/r ps = (1/t p „)<? occuring in the denominator of the 
w-integral in the expression for (see Eq.(4) in Ref. 8) is comparable to or larger than the other scattering rates 
due to phonon-defect scattering, we expect that K p h rapidly decreases for increasing field for temperatures smaller or 
close to T ~ (1/2) A(T, H) . In the measurements of k one observes a fast drop of k in small fields at temperatures 
around T c /6 4 which can be attributed to the fast drop of n p h due to scattering by quasiparticles in the 7r4^and vortex 
lattice. 

It is interesting to compare the temperature dependence of the phonon lifetime r p with the lifetime r e of a quasi- 
particle due to scattering by acoustic phonons. This scattering rate has been derived for the BCS state in Ref. 19 
(see Eq.(3.9)). We extend this expression for the scattering rate l/r es in anology to our procedure for 1/t p from the 
BCS to the vortex state by replacing the density of states by the DOS A(Q) in Eq.(l) and the coherence factor terms 
by B(tt) in Eq.(2). Then we obtain: 

l/r„(w) - C[/(-/k>)]~ lj4_1 ( fi ) / dLJ w^{A(n)A(fl-n ) - B(fl) B(fl - fl ) 

Jo 

xf[-0(u-e;o)][l + b(0u> o )] 

+ [A(fl) A(fl + fl ) - B{Q) B(fl + fi )]/[-/3(w + w )]&(/3w )} . (15) 

Here C is a constant arising from the electron-phonon matrix elements, (3 = 1/T, b(x) = (e x — 1) _1 , and ojq is 
the phonon frequency. We use this expression now in the calculation of the electronic thermal conductivity n e . We 
mentioned above that the integrand of the w-integral in the expression for n e given in Ref. 7 contains the term 
Ime = 7 + 7.4 in the denominator which corresponds to the sum of relaxation rates of a quasiparticle due to impurity 
and Andreev scattering. The effect of thermal scattering by phonons can now be taken into account by adding 
l/r es to Imeo • The temperature dependence of n e is governed by the factor (w/T) 2 sech 2 (aj/2T) in the integrand of 
the w-integral which yields the main contribution of the integrand from the vicinity of x = oj/T at about x m = 2 . 
Therefore we can estimate the temperature dependence of 1 /r es by carrying out the following variable transformations 
in Eq.(15): x = uj/T , t = ujq/T , £1 = xy, and £Iq = ty with y = T/A(T,H) . This yields l/r es as a function of 
il = xy and x. The normal state relaxation rate l/r e „ is obtained from Eq.(15) by setting A = 1 and B = 0. In 
Fig. 5 we have plotted r en /T es as a function of y for constant x = x m = 2 and A = 0.2, 0.3, and 0.6 (H/H C 2= 0.78, 
0.62, and 0.29, for S = 0.1 , see Fig. 1). One sees that for increasing A, or decreasing field, this function tends to a 
step function where the step occurs at y = 1/2 , that is, at a quasiparticle energy Q = oj/A = 1 . This result agrees 
with the previous result for the BCS state. 19 One recognizes that the variations of T e „/r es as a function of y are much 
smaller than those of the phonon lifetime ratio T pn /r ps shown in Fig. 4 where, for increasing A , a step occurs at 
about y = CI/ A = 1/2, i.e., at a phonon energy Oo = u>o/A = 2 . The contribution of l/r es to the relaxation rate Imeo 
occuring in the electronic thermal conductivity K e increases with temperature like T 3 because 1 /r en is proportional to 
T 3 . In the presence of moderate impurity scattering this becomes comparable with 7 only at temperatures such that 
y = T/A(T,H) > 1 where, according to Fig. 5, the ratio T en /r es is nearly one. The expression for l/r es in Eq.(15) 
has been calculated in analogy to Ref. 19 for quasiparticle scattering by acoustic phonons. Instead of this spectrum, 
ClOq , one should employ in the expression for l/r es the transport Eliashberg functions al r F nn {uLi ) , where n — a, n , 
which have been calculated for MgB2 . 21 Since the coupling constant X aa = 0.81 > = 0.41 , we conclude that 
for the cr-band the field dependence of l/r es should be taken into account while, for the 7r-band, l/r es can safely be 
approximated by 1 /r en . 
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III. CONCLUSIONS 



We have calculated the density of states A(u) , the electronic thermal conductivity n e , and the scattering rates 
1/t p (uj) of a phonon and 1/T e (w) of a quasiparticle for a single isotropic gap over the entire field range of the mixed 
state. That the BPT 9-11 and P 12 methods for these calculations arc valid approximation schemes has been proved 
by comparison of the results for A with the solutions of the Eilenberger equations for Abrikosov's vortex lattice order 
parameter. 18 By solving the gap equation including impurity self energy and vertex corrections (see Fig. 1), we 
have shown that the relation between the central parameter of the theory, A = A A/t> , and H/H C 2 can be very well 
approximated by the Ginzburg-Landau (GL) relation [see Eq.(5)] over the entire field range if one uses an appropriate 
value of the Abrikosov parameter Pa ■ This result is connected with the fact that the Abrikosov order parameter is 
a solution of the linearized GL equations which allows extension of the BPT method over the entire region of linear 
magnetization. 22 

The second main result of the present paper is that the vertex corrections due to impurity scattering 10,12 substan- 
tially reduce the dependence of the functions r y s (H)/j n , K es (H)/n en , and T pn /T ps (H) on the amount of impurity 
scattering. The impurity scattering rate 8 = T/Aq = (tt/2)(^o/£) is renormalizcd here in the Born approximation 
r — > 7 = r A(tS) where A(u>) is calculated self-consistently from Eq.(l). The vertex corrections yield a renormalization 
of the gap A r due to the function D [see Eq.(7)]. As can be seen from Fig. 2, the functions ') s {H)/ r ) n for 8 =0.1, 0.5, 
and 1.0 are not very different from each other. The results are however quite different if we set the renormalization 
function D equal to one. Then the slope of j s (H)/j n versus H/H c2 at H = rises rapidly with 8 and becomes very 
large for 8 = 1 . It should be noticed that all curves in Fig. 2 tend to zero in the zero field limit as they should to be 
in accordance with Anderson's theorem for "dirty" superconductors. Our result that the effect of impurity scattering 
on the field dependence of "/ s (H) /"/ n for a single isotropic gap is rather small and does not yield a rapid rise for small 
fields, gives support for the two-gap model in interpreting the measurements oi ^f s (H)/j n on single crystals of MgB2 . 3 
According to this model, the first steep increase of 7 S (H ) /"f n , which is nearly the same for applied fields parallel and 
perpendicular to the c-axis, arises from the behavoir of the small gap associated with the three dimensional 7r-band. 
This in effect closes for increasing field due to the overlap of the vortex cores with large radius £ w <~ 500 A 6 at a 
relatively small "virtual" upper critical field H^ 2 ~ 3 - 4 kOe. The remaining superconducting contribution above H^ 2 

(c) 

arises from the large gap associated with the two dimensional cr-band which closes at an upper critical field H c2 ~ 
32 kOe and ^160 kOe. The field dependence of the 7r-band and the cr-band contributions to j s (H)/^f n is 

presumably given by the almost linear curves in Fig. 2 for 8 = 0, 8 and 8 = 0.5 with H c2 equal to H^ 2 and or 
H^ 2 b ^ , respectively. 

The results for n es (H) / n en shown in Fig. 3 depend on the impurity scattering rate 8 primarily because the effect of 
Andreev scattering on the upward curvature becomes relatively smaller as 8 increases. Although the slope at H = 
is relatively small even for 8 = 1 , it becomes very large at 8 = 1 if the renormalization of the gap A r — A/D due 
to the vertex correction D in n e [see Eq.(13)] is neglected. This new result has the consequence that the observed 
steep rise of K e (H) for small fields parallel or perpendicular to the c-axis in single-crystalline MgB2 4 can again be 
explained only in terms of the two-gap model where the smaller gap in the three dimensional 7r-band closes in effect 
at a "virtual" upper critical field Hf 2 ^3-4 kOe due to the overlap of the vortex cores with large radius £ n . 6 At 
higher fields n e {H) tends to saturate until, in the vicinity of ~32 kOe, the contribution to K e from the larger 
gap associated with the cylindrical cr-band exhibits an upward curvature. The contributions to n e arising from the ir- 
and cr-bands are presumably given by the curves in Fig. 3 for 8 = 0.8 and 8 — 0.5 with the corresponding i/ c2 equal 
to Hf 2 and where the curve for (5 = 0.5 yields the measured upward curvature near H C 2 . 

The rapid decrease of k observed at low fields, which is attributed to the decrease of the phononic thermal con- 
ductivity K p h , 4 can be explained by the field and relevant temperature dependence of the ratio of phonon lifetimes 
T pn I 't P s(H) — g(H) shown in Fig. 4. One sees that g(H) tends rapidly to the constant one for increasing field and 
values of y < 1/2 corresponding to temperatures T < (1/2)A(T, H) . If we assume that A is the smaller gap associated 
with the 7r-band, then the T satisfying this condition are much smaller than T c . However, at the lowest temperatures, 
this effect on n p h is strongly reduced because the scattering rate 1 / r ps occuring in the denominator of the w-integral 
for Kph is multiplied by T/T c in comparison to the scattering rate due to sample boundary scattering (see Eq.(4) of 
Ref. 8). 

We have calculated n e (H) in the limit T — > , or uj — > . For finite T one has to solve the transcendental equation 
determining the pole £o of the BPT-Green's function 7 for each finite value of u> occuring in the w-integral for K e . It 
should be pointed out that this calculation is avoided in the expression for n e which was derived from linear response 
theory 13,14 For higher temperatures one should add to the sum of scattering rates 7 and 7^4 due to impurity and 
Andreev scattering [see Eq.(13)], the scattering rate l/r es (7J) due to scattering of quasiparticles by phonons. The 
field and relevant temperature dependence of the ratio T en /T es (H) is shown in Fig. 5 as a function of y = T/A(T, H) . 
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One sees that this ratio is nearly equal to one for values of y > 1 corresponding to temperatures T > A(T,H) . For 
the small gap associated with the 7r-band this condition is probably satisfied at those temperatures where 1 /r e „ ~ T 3 
becomes comparable to the intraband impurity scattering 7 for the 7r-band. However, for the large gap associated 
with the CT-band, this condition might not be satisfied. We conclude from these estimates that, for the 7r-band, l/r es 
can be approximated by l/r e „ while, for the er-band, the field dependence of l/T es (H) should be taken into account. 
These considerations are supported by detailed calculations which yield a larger intraband impurity scattering rate 7 
for the 7r-band than for the cr-band, 20 and a larger scattering rate 1/T en for the cr-band than for the 7r-band because 
the intraband electron-phonon coupling constant satisfies \ aa > . 21 

Our assumption that we can simply add our results obtained for two different isotropic gaps and different upper 
critical fields in the ir- and c-bands is a crude approximation in view of the actual situation in MgB2 . It is true that the 
interaction due to interband impurity scattering between the a- and 7r-bands is small. 20 However, we have neglected 
the pairing interaction between the two bands which presumably leads to Cooper pair tunneling from the cr-band with 
strong pairing to the 7r-band where it leads to giant vortex cores and where superconductivity is maintained above 
the virtual H% up to . 5 Nevertheless we believe that our clear-cut results for the mixed state of a single isotropic 
gap are helpful in analyzing the complicated behavior of different physical quantities in MgB2 . 

Another approximation is that we have neglected higher harmonic oscillator components N of Abrivosov's vortex 
lattice order parameter which lead to a distortion of the lattice at lower fields. 18 For an s-wave superconductor the 
BPT method breaks down at low fields where properties are determined by the states bound to the vortex cores. 
However, STS measurements in MgB 2 show an absence of localized states in the cores. 6 
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Fig. 1 
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FIG. 1. Solutions of gap equation for H/H C 2 vs A = AA/u (A = (2eH)^ 1 ^ 2 ), for impurity scattering rate 8 = 0.1, 0.5, and 
1.0, for 9 = 7r/2 (solid curves). The fits are obtained from the analytical expression with the Abrikosov parameter Pa = 1.1, 
1.4, and 1.7 [see Eq.(5)]. The dashed curves show the #-angle averages for 8 = 0.1, 0.5, and 0.8 with fit parameter values /3a = 
2.3, 2.8, and 3.2. 
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Fig. 2 
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FIG. 2. Ratio of specific heat coefficients, 7s/7n , vs H/H C 2 for reduced impurity scattering rates S = T/Ao = 0.1, 0.5, and 
1.0 (solid curves) for 6 = 7r/2 (cylindrical Fermi surface(FS)). The dashed curves show the #-angle averages (spherical FS) with 
8 and Pa as in Fig. 1. 
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Fig. 3 




FIG. 3. Electronic thermal conductivity ratio n es /K en vs H/H C 2 for 8 = -k/2 (solid curves) and for the #-angle avera 
(dashed curves). Notation as in Fig. 2. 
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Fig. 4 
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FIG. 4. Ratio r pn /r ps of phonon lifetimes vs y — T/A(T,H) corresponding to dominant phonon energy ojo = 4Aj/ = 4T 
which determines the phononic thermal conductivity k p h ■ The constant fields are H/H C 2 = 0.29, 0.62, 0.78, or, A = 0.6, 0.3, 
0.2 (from bottom to top) and S = 0.1 . 
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Fig. 5 




FIG. 5. Ratio r en /r es of quasiparticle lifetimes due to phonon scattering vs y = T/A(T,H) corresponding to dominant 
quasiparticle energy to = 2Ay = 2T which determines the electronic thermal conductivity K e . The constant fields are H/H C 2 = 
0.29, 0.62, 0.78, or, A = 0.6, 0.3, 0.2 (from bottom to top), for 8 = 0.1 . 
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